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Linear programming techniques are becoming of greater 
importance because the use of computerization has increased 
the fields for applications for linear programs. The primal- 
dual algorithm, in which the constraints are added one at a 
time, is investigated as a possible faster solution method. 

A computer program was developed to compare this method with 
the standard primal-dual algorithm using the full set of 
constraints at one time. Several random problems were solved 
using these two methods, and the results indicated a signif- 
icant improvement in the solution time by the use of adding 
the constraints one at a time. 
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I 




INTRODUCTION 



I . 

New linear programming algorithms have been developed to 
reduce the computational time in solving linear programs . The 
purpose of this thesis is to investigate the merits of one 
such new algorithm. This method consists of introducing the 
constraint equations one at a time. After each constraint is 
added, the "smaller", or submatrix, problem is solved using 
the primal-dual algorithm. This continues until all con- 
straints have been added and a solution is obtained. 

The rationale for this approach is that small matrices 
are used in the initial stages of solving the linear program; 
the size of the matrices increases only when additional con- 
straints are introduced. If the number of iterations used in 
this method is not significantly different from the number of 
iterations used with the full matrices, the manipulation of 
the smaller matrices in the initial stages will reduce the 
solution time. 
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II. NOTATION 



m 

n 

A 




P . 
ai 



X . 

3 

C 



B 



s . 
3 

/N 

S . 

3 

d 



number of constraint equations, 
number of legitimate variables. 

mxn matric of coefficients of the constraint equations 
with elements j / i = l/.../m, j = l,...,n . 
mxm matrix of the basis vectors, 
inverse of the basis. 

mxl column vector which is the j column of A. 




t 3 lf.../n . 

th 

mxl column vector associated with the i^ artificial 
variable, i = l,...,m . 

mxl column vector with elements x. where X. = P ^P . . 

1 ] 3 3 

nxl column vector with elements c . which are the costs 
of the legitimate variables. 

mxl column vector with elements b^ which are the right- 
hand sides of the constraint equations, 
dual slack variables. 



i . e . , P . = 
3 






dual slack variables after a dual iteration. 



mxl column vector whose elements are the coefficients 
of the basis variables of the added constraint 
equation. 
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Notation used with tableau; 





-4/-1 



© 




the vector to be introduced into the basis, e.g. , 
in this tableau Pq will be introduced, 
pivot element for a primal iteration, i.e., the 
6 criterion is 6 = min(x. /x. .) such that x. . > 0. 



pivot element for a dual iteration, i.e., the 

6-criterion is 6 = min (-s ./z --c . ) such that 

j 3 13 



z .-c . < 0 . 
D 3 
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III. FORMULATION OF THE PROBLEM AND SOLUTION PROCEDURE 
The general linear programming problem is to maximize 

n 

z = y c .X . 

jil ^ 

subject to 



n 



y a^jXj ^ ~ l/.../m , 



j = l 



and 



Xj ^ 0, j = l,...,n 



( 1 ) 



The modified primal uses an additional constraint 



X, 



n 

j-1 J 



“ / 



where the cost of Xq is zero and b^ is arbitrarily large, so 
that for Xq > 0 the constraint adds no additional restriction 
on (1) . 

The modified primal is written to maximize 

n 

c .X 

D 



n 

= y c . X . 
jii 



subject to 



n 



J=1 ■> 



n 

y a..x.=b.,i=l,...,m, 
.t. ID D 1 



( 2 ) 



j = l 



and 



Xj ^ 0 , j = i , . . . ,n . 

From (2) we can write the modified dual with slack var- 
iables, Sj, j = 0,l,...,n, added. That is to minimize 

m 

w,,b,. + y w.b. 

0 0 . ^ , 11 
1=1 
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subject to 



w 



0 



= 



0, 






m 

I 

i=l 



w . a . . 
1 1 ] 



- s . = c 



j = 1,. 



.n, 



(3) 



and 

unrestricted for i = 0,1,..., m. 

The starting feasible solution to the primal dual al- 
gorithm is = 0 , i = 1 , . . . ,m, and Wq = max (c . , 0) . 

^ j ^ 

For an optimal solution the complementary slackness con- 
dition must hold. That is 



n 

3=1 -* 



Adding artificial variables, x ., i = 0,1,..., m, to (2) 



with the cost of the artificial variables set to -1 and the 



cost of the legitimate variables set to zero, the extended 



primal can be written as 
maximize 

m 

- ’‘ao - .1 ’‘ai 
1=1 



subject to 

n 

■'o + .1 ==30 = ‘■o ' 

3=1 -J 



n 

y a..x. + X . = b., i = l,...,m, 
ID 3 1 



and 



X. > 0, j = l,...,n, and x . > 0, i = 0,...,m . 
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Solving the extended primal is similar to using a Phase 
I Revised Simplex method. However, in the primal-dual 

algorithm, when Phase I ends the linear program is solved 
because complementary slackness is maintained throughout the 
solution procedure. 

If a new constraint is added to the tableau, complemen- 
tary slackness is maintained without changing the dual slack 
variables. This can be shown as follows: 

Assume we have a feasible solution to the problem 
with k constraint equations. This means that s^ = 0 

for all j such that P. € P and x. = 0 for all j such 

3 3 

that Pj ^ P. These conditions imply that complemen- 
tary slackness is maintained, and that the modified 

dual also has a feasible solution. 

s t 

Now we add the k + 1 constraint which introduces 
a new dual variable, w, . , , but no new dual slack 
variables, s ^ , j = 0,1,. ../n. We need a feasible 
solution to the modified dual for the enlarged 
system. Observe that we have a feasible solution if 
we set = 0 since then the s ^ , j = 0,1,..., n, 

remain unchanged. In particular, s^ = 0 for all j 
such that Pj e P, that is, for the legitimate var- 
iables. Also, Xj = 0 for all j such that P^ ^ P, 
which implies that we have maintained complementary 
slackness . 

It is worth noting that the z ^-c^ , j = 0,1,..., n, must be 
recalculated since the new constraint which is added to the 
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extended primal starts with its artificial variable, x vj.i ' 

3. f Kt i. 

in the basis with its cost set at -1. 

P 0 0 

The new basis is P = , where P , , = 

_m ^ f 

d^ 1 1 

is the artificial vector associated with x , , , . Now we can 

a,k+l 

solve the new k + 1 system using the primal-dual algorithm 
since complementary slackness has been maintained. 

An optimal solution exists if and only if the following 
criteria are satisfied: 

1. z.-c. > 0 for j = 0,1,. ..,n, 

D D - 

2. z^-c^ = 0, and 

3. Xq > 0. 

The solution procedure is as follows: 

The first tableau is set up using the first two 

constraints of the extended primal and the starting 

solution to the modified dual, which implies that 

at least one s. = 0. 

3 

Step 1. Is there a j, say jQ, such that s^^ = 0 and 



jo 



"j0-=j0 •= 

a. Yes, Go to 2. 

b. No. Go to 3. 

Step 2. Introduce P^q into the basis using the minim\im 

0-criterion and a primal iteration. Since the extended 

primal is bounded a pivot will always exist. Note that the 

Sj remain unchanged for all j. Go to 1 . 

Step 3. Is z.-c. < 0 for some i? 

3 3 

a. Yes. Go to 4. 

b. No. Go to 5. 
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Step 4. 


Use the minimum 0-criterion. Is 0 bounded? 


a. 


Yes. Perform a dual iteration to compute a new 







set of s.'s, say s.'s. Go to 1. 

3^3 

b. No. The linear program has no feasible solution. 



Stop . 




Step 5. 


Is Zg-Og < 0? 


a. 

Stop. 


Yes. The linear program has no feasible solution. 


b. 


No. Go to 6. 


Step 6. 


Have all of the constraints been added? 


a. 


Yes. Go to 9. 


b . 


No. Go to 7. 


Step 7. 


s t 

Introduce the next restraint, say the k + 1 



Place the artificial vector P , , , in the basis. Compute 

a , k+ 1 ■ 



^B,k+1* 


=B,k+l i “• 


a. 


Yes. Go to 8. 


b. 


No. Multiply all coefficients of the new constraint. 



except for the artificial variable, by -1. This assures that 
Xg ^ 0 and the artificial variable is non-negative. Go 

to 8 . 

Step 8. For the system with k + 1 restraints, compute the 

new values of z.-c. for j = 0,1,..., n. Go to 1. 

3 3 



Step 9. 


Is Xq = 0? 


a. 


Yes. The linear program is unbounded. Stop. 


b. 


No. An optimal solution has been found. Stop. 
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IV. SAMPLE PROBLEM 
Consider the following example: 
maximi ze 



z = 2Xj^ + 4x^ 

subject to 

3x^ + 4 X 2 + 6 X 2 + = 24 , 

4Xj^ + 3 X 2 + 12x^ + x^ = 24 , 

Xi + X 2 + 4x^ = 8 , 

and 

Xj ^ 0, j = 1,...,5. 

Then the extended primal is 
maximize 



-x»-x, - x»-x, 
aO al a2 a3 



subject to 



+ x^ 


+ 


X2 


+ 


x^ + X 


3 Xi 




4x2 




6X2 + X 


4 xi 


+ 


3x2 




12 x 2 


^1 




^2 




4X2 



+ X 



al 



+ X. 



+ X 



a 2 



= bQ, 

= 24, 
= 24, 



+ X , = 8 , 

a3 ' 



X. > 0 for j = 0,,..,5, and x . > 0 for i = 0,...,3 . 

J — J r r r ai — 



The dual slack variables are Sq = max(c.,0) = 4 with j = 3. 

j ^ 



Then s„ 

13 

that Sj^ 



x,.b_ = 4b», and s. = s„ - c. for j = 1 
0 0 0 2 0 •) 

2, S 2 = 4, = 0; = 4, and s^ = 4. 



5 , so 
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The starting tableau, using the first original constraint, 



is 




24/6 = 4 



From step 1, we see that = 0 and z^-c^ < 0. Using 

the minimiam 6-criterion (as discussed in section III) in step 

2, we introduce into the basis and remove P , from the 
' 3 al 

basis . 



Since there is no j,. for which s.,, = 0 and z.^-c.- < 

-^0 3O 3O 3O 

but z.-c. < 0 for several values of i, we arrive at step 4 

3 3 j/ M 

Using the minimum 0-criterion a new set of Sj'®/ called s^ 
are calculated. 



0 , 



's 



p 


B 


^0 


^1 


^2 


^3 


P4 


^5 


^aO 


^al 


^aO 


1 

0 

XJ 


© 


1/2 


1/3 


0 


5/6 


1 


1 


-1/6 


^3 


4 


0 


1/2 


2/3 


1 


1/6 


0 


0 


1/6 


z.-c. 
3 3 


-bo+4 


ED- 


■1/2 


-1/3 


0 


-5/6 


-1 


— 


— 


s . 
3 




4 


2 


4 


0 


4 


4 


— 


— 


/\ 

s . 
3 


16 


0 


0 


8/3 


0 


2/3 


0 


— 


— 



Now Sq= 0 and z^-Cq < 0 so, from step 2, 
Pq into the basis and remove P^q from the bas 



e = 



bo-4 



0= -4/-I = 4 



we introduce 
is . 
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p 


B 


^0 


^1 


^2 


^3 


^4 


^5 


^aO 


^al 


^0 


1 

o 


1 


1/2 


1/3 


0 


5/6 


1 


1 


-1/6 


^3 


4 


0 


1/2 


2/3 


1 


1/6 


0 


0 


1/6 


z . -c . 
3 3 


0 


0 


0 


0 


0 


0 


0 


— 


— 


s . 
3 


16 


0 


0 


8/3 


0 


2/3 


0 


— 


— 



From the above tableau we trace through steps lb, 3b, 

5b, 6b, and arrive at step 7. Note that we have obtained an 

optimal solution to the subproblem with one constraint. In 

step 7 we introduce the second constraint. Since Pq and 

-T 

are basis vectors, d = (a 2 Q,Si 22 ^ ~ (0^2); the new basis 
consists of Pq, P^ and P^ 2 * With this basis we find that 
x„_ = -24 < 0 so that step 7b must be used. The- second re- 
straint is replaced by 



-4^ - 3„-12--X|- + x„ = -24 , 
xl x2 x3 5 a2 ' 

which is used throughout the remainder of the solution proce- 
dure. Note that now d"^ = (0,-12) and Xg 2 = 24 > 0. New 
values of ^j“^j 3 = 0,1,..., n are computed (step 8), and 

the new tableau is 




1/2 



= 8 
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since = 0 and < 0 (step 1) we go to step 2. 

Using a primal iteration, we introduce into the basis and 



eliminate from the basis. 



p 


B 


^0 


^1 


^2 


^3 


^4 


^5 


^aO 


^al 


^2 


^0 


00 

1 

o 


1 


0 


-1/3 


-1 


2/3 


1 


1 


-1/3 


0 


^1 


8 


0 


1 


4/3 


2 


1/3 


0 


0 


1/3 


0 


^a2 


8 


0 


0 


7/3 


-4 


0 


-1 


0 


4/3 


1 ( 


z.-c. 
3 3 


-8 


0 


0 


-7/3 


4 


-4/3 


1 


— 


— 


( 


s . 
3 


16 


0 


0 


8/3 


0 


2/3 


0 


— 


— 


— 


/S 

s . 


12 


0 


0 


3/2 


2 


0 


1/2 


_ _ _ 


— 


... 


3 


























= 6 



-2/3 _ / 



Now z.-c. > 0 for all j for which s. =0. From steps lb, 
3 3- 3 

3a, and 4a, a new set of s^'s are computed. Then s^ becomes 
zero and z^-c^ < 0 so , from step 2, P^ enters the basis and 
P ^2 removed from the basis. 




Optimality has now been obtained with the second restraint 
added. Following steps lb, 3b, 5b, 6b, and 7, we introduce 
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the third and final restraint. Since the basis vectors were 

-T 

Pq, P^, and P^, d = (^ 39 ' ^ 31 ' ^ 34 ) ~ (0,1,0). The new basis 
vectors are Pq, P^, P^, We find that = 2 > 0 so we 

go to step 8 and recompute z^-c^ for j = 0 , 1 , ...,n . 

The next sequence of steps is lb, 3a, and 4a, which 




Steps la and 2 bring P^ into the basis with the elimina- 
tion of P o from the basis. The new tableau is: 
a3 
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This tableau is the final tableau since the sequence of 
steps lb, 3b, 5b, 6a, and 9b inform us that we have found an 
optimal solution to the original linear program. Note that 
the complementary slackness condition has been maintained, 

n 

that is, s^Xp, + 'l s.x. = 0. Therefore, the optimal solution 

u u j^l 3 3 



is 



X 



0 




14 > 0, 



''l = ''2 = ""S = 
Xj = 2, and 



X4 — 12, 



with the optimal cost z = Sg = 8 . 
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V. PROGRAMMING TECHNIQUE 



The linear programming technique described in this thesis 
was programmed in FORTRAN IV for use on the IBM 360/67 com- 
puter. One subroutine is used for the primal simplex itera- 
tion and another is used for the dual iteration. The final 
subroutine is used for the addition of constraints. The main 
(driving) routine is used to solve both the full linear pro- 
gram and the linear program using addition of constraints as 
described in this thesis. Using the main routine to solve 
both problems eliminates any time differences due to differ- 
ences in programming techniques. Both of the solution proce- 
dures were timed*, and the number of iterations of each were 
counted. Read and print times were not included in the 
timing . 



*The timing routine was developed by Lt. E.A. Singer, a 
student at the Naval Postgraduate School. 
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VI. EFFICIENCY OF THE ALGORITHM 



To eliminate the considerable time and effort required 
to input data by hand, a subroutine was designed which gen- 
erates random problems of a large size. This routine uses a 
random number generator to produce the elements of the A, B, 
and C matrices. The following criteria were used in order 
to insure the existence of a bounded feasible solution: 

maximize 

n 

z = y c . X . 

jil 

subject to I 

n 

y a..x. - X . = b., i = l,...,m, 
ij D SI i' 

and 

c . < 0 , X . > 0 , X . > 0 , a.. > 0, b. >0 
3 - J - SI- 13- 1- 

f o r r 1 , . . . , m , ^ 1 , . . . , n , 

th, 

where x . is the slack variable for the i constraint. 

SI 

For this investigation the subroutine generated problems 
having 70 variables (including 20 slack variables) , 20 con- 
straint equations, and 50 cost coefficients using the follow- 
ing uniform distributions: 

a^^j, uniform (0,1); 

bj^ , uniform (0,5) ; 

Cj, uniform (-1,0). 

A total of 40 problems were solved using the random 
problem generator described above. The execution times for 
these problems are given in Table I at the end of this section. 
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The comparison of solution times shows that ail 40 prob- 
lems ran faster using the method described in this thesis, 
than with the standard primal-dual algorithm. The time dif- 
ferences range from 15,17 seconds to 48.18 seconds with an 
average time difference of 27»53 seconds. 
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TABLE I 



Prob . 
No. 

1 


Full Array 


Addition of 
Constraints 


X . - y . 
1 -^1 


Time (x . ) 
(sec. ) ^ 


Iter. 


Time (y . ) 
(sec.)^ 


Iter 


1 


87.82 


23 


56.50 


21 


31. 32 


2 


87. 84 


23 


56.91 


23 


30.93 


3 


80.22 


21 


56.53 


21 


23.69 


4 


87.93 


23 


56.70 


22 


31.23 


5 


84.06 


22 


58.61 


23 


25.45 


6 


91.63 


24 


56.62 


22 


35.01 


7 


95.45 


25 


58.65 


23 


36.80 


8 


80.31 


21 


56.52 


21 


23.79 


9 


99.60 


26 


65.67 


32 


33.93 


10 


80.16 


21 


56.51 


21 


23.65 


11 


80.18 


21 


56.50 


21 


23.68 


12 


91.65 


24 


65.12 


29 


26.53 


13 


80.19 


21 


57.99 


22 


22.20 


14 


91.74 


24 


64.35 


30 


27.39 


15 


80.20 


21 


56.48 


21 


23. 72 


16 


84.02 


22 


57.10 


22 


26.92 


17 


95.52 


25 


59.36 


24 


36.16 


18 


83.95 


22 


58.71 


22 


25.24 


19 


84.14 


22 


65.47 


25 


18.67 


20 


99.25 


26 


70.03 


28 


29.22 


21 


107.08 


28 


58.90 


25 


48.18 


22 


87.85 


23 


56.51 


21 


31.34 


23 


84.01 


22 


61.39 


24 


22.62 


24 


80.16 


21 


56.52 


21 


23.64 


25 


80.17 


21 


56.49 


21 


23.68 


26 


80.16 


21 


56.50 


21 


23.66 


27 


80.19 


21 


56.52 


21 


23.67 


28 


84.21 


22 


58.84 


23 


25.37 


29 


103.27 


27 


57.85 


24 


45.42 


30 


87.83 


23 


59.47 


23 


28.36 


31 


84.07 


22 


59.33 


22 


24.74 


32 


88.19 


23 


56.92 


23 


31.27 


33 


80.28 


21 


56.56 


21 


23.72 


34 


87.96 


23 


56.56 


21 


31.40 


35 


84.11 


22 


57.09 


23 


27.02 


36 


84.05 


22 


56.56 


21 


27.49 


37 


80.20 


21 


56.51 


21 


23.69 


38 


80.18 


21 


56.54 


21 


23.64 


39 


84.03 


22 


68.86 


26 


15.17 


40 


80.29 


21 


56.54 


21 


23.75 




Total 


3454.01 




2352.77 




1101.24 


Average 


86.35 




58.82 




27.53 
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VII. SUMMARY AND CONCLUSIONS 



A modification of the primal-dual algorithm has been 
presented. This modification differs from the standard 
primal-dual algorithm in that the constraint equations are 
introduced one at a time, and each subproblem is solved be- 
fore the next constraint is added. 

This algorithm was programmed in FORTRAN IV for the IBM 
360/67. The program was designed so that any given linear 
program is solved first by the standard primal-dual algor- 
ithm, and then is resolved using the modified primal-dual 
procedure. Further, the same subroutines are used for both 
methods in order to eliminate timing bias due to coding 
differences. In fact the modified procedure contains steps 
which are not included in the timing of the standard routine. 

In all cases the modified procedure was faster than the 
standard procedure. Of 40 test problems the standard method 
averaged 86.35 seconds per problem, whereas the modified 
method averaged 58.82 seconds per problem. One should not 
judge the actual running time of the test problems since no 
attempt was made to improve the efficiency of the computer 
program on an absolute basis; only the relative speeds of 
the two methods is of importance. 



23 



BIBLIOGRAPHY 



1 . 



Hadley, G. , Linear Programming . Reading, Mass: 
Wesley Publishing Co., Inc., 1962. 



Addison- 



24 



APPENDIX A 



FLOW DIAGRAMS OF THE COMPUTER PROGRAM 
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1 . MIIIN PROGRAM 





/ V/r wir^' 

TPei.£/^U 



77)3lP/?i/ J 



call 

, SrMPLX ( j,y 



yVo-r£- 




//0 J 




V‘j. 



^(-%1 \ -s';. 



P^/Z gj -c/<o 

1 S P^/lj/^A£J> 



^ piie P 7 



CALL 



\ 



(/) Iu/vS= O) f^B^^X8l^£ ioV 

£KI^T^ (c<f>A/TI/^H£ ) > ' f 

(Z) JUhy 8 ^ / ) Uh/8^Uf/PiX> StfLUTT^^ ^ T 

£')fl£T‘^ J 3UT AQpP£ C^M/^PPPIPTSI ' ^ ^ ^ 

ry 8£ ADD££>» 1 M P L X ( j ^)) 

C3) IUf/6- 2 '^ ' Uh/a^UA/P£i> ^^XM7I^A/\ j 

pyi^r^ aa/X> A 7.1. c^A>3'rPA£^rJ' 

h/py’£ g£E// A£>£iE £> . 
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SIMPLEX ITERATION SUBROUTINE 




28 



DUAL ITERATICN SUBROUTINE 




-4. ADDITION OF NEW RESTRAINT SUBROUTINE 




A-a/£) 



7 



r ■•” 
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APPENDIX B 



FORTRAN LISTING OF THE COMPUTER PROGRAM 
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0001 

OOC2 

0003 

000^ 

cocs 

0006 

0007 

OOCB 

000^ 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

OCIO 

0020 

00 21 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0020 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

003Q 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0C51 

C052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 



ooo noo ooo nr>r> 



IV G LEVEL 1, MOO 1 



MAIN 



DAT'= 



6815^^ 



05/22/31 



■ - - ‘ 0000040 

0000045 



COMMON/ I NT/ I UNR, I M , M , N , I R A S » I T ER 
COMM0N/FLQAT/P,H,C, RT,0NFfTEM,PlNV,7CR, 7 C , W S B , WS , T , X , ORAR , FP 



NOTE - t^S(J) IN THIS PROGRAM ANO S(J) 

EP = lO.OE-6 
XF = URN( 0) 

IPRQB = 0 

6 IF( IPROB.GE. 50» GO TO 50C 
IPROB = IPR('B +1 
call PSEUf'tn 

M = M 4- 1 
N = N + 1 



IN THESIS ARE EQUIVALENT 

0000050 

0000055 

0000060 

0000070 

0000080 

0^00080 

0000000 

0000100 



KO=l IS FOR USE OF COMPLETE TABLEAU 

KO = 1 
IM = M 

INITIALIZE TABLEAU 

5 on 1 I= 2 fM 
BT( I , 1 ) = O.C 
BT( I , 2 ) = B( I-l ) 

1 P( I, 1 ) = 0.0 
ITER = 1 

on 2 J= 1 ,N 

2 P( If J) = 1.0 
0 T( 1 , 1 ) = 1.0 
3 T( 1 , 2 ) = O.C 
00 3 1 = 1 , M 
IBAS( I ) = N+I 

3 ONE( I) = -l.C 

on 10 1 = 1 ,^ 

DO 10 J= 1 ,M 
IF( I .ED. J) GO TO 9 
PINV( I , J) = 0.0 
GO TO 10 

9 P INV( I , J) = 1.0 

10 CONTINUE 

CALL TI*^EIT(C,CL) 

7 CB( 1 ) = - 1.0 
ZC 3 ( 2 ) = 0.0 
00 11 1 = 2 , IM 

11 ZrB( 2 ) = ZC 4 ( 2 ) -BT( I , 2 ) 

ZC( 1 ) = - 1.0 

on 20 J = 2 ,N 
ZC(J) = -l.C 
DO 20 I = 2 ,TM 
20 ZC( J) = ZC ( J » -P(! , J) 

INO = 0 
WS( 1 ) = 0.0 
DC 30 J= 2 ,N 

IF(C( Jl .LE.WS( 1 ) )G 0 TO 30 
WS( 1 ) = C( J) 

INO = J 
30 CONTINIJF 

wSB( 1 ) = VIS( 1 ) 

WSB( 2 ) = O.C 
DO 40 J= 2 ,N 
40 WS( J) = WS( 1 ) -C( J) 

OFTERMINE VECTCR TO INTRODUCE 

IF( INO.EO.O.OGO TO 50 
IF( ZC( IND» .GF.O.OIGQ TO 70 
CALL SIMPLX(IND) 

GO TO 60 

50 CALL SIMPLXI 1 ) 

60 IF( IUNB- 1 ) 7 C, 130,190 
70 TFMP = 9 Q 9999 . 



0000260 

0000265 



0000770 

0000280 

0000790 

0000300 

000^305 

0000310 

0000370 

0000375 

0000330 

0000340 

0000380 

0000360 

0000370 

0000380 

00003Q0 

0000400 

0000410 

0000470 

0000430 

0000435 

0000440 

0000450 

0000460 

0000470 

0000480 

00004QO 

0000500 

0000510 

0000570 

0000530 

0000540 

0000550 

000056C 

0000570 

0000580 

00005Q0 

0000600 

0000610 

0000620 

0000630 



000064^ 

0000650 

0000660 

0000670 

0000680 

0000690 

0000^00 
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FORTRAN IV G LEVEL 1, MOO 1 



MAIN 



DATE 



6B159 



05/22/31 



£ FIND THETA FOR DUAL = MI N ( -S ( J ) / ZC (J )) OVER J FOR ZCUXO 
C 



0061 






IND = 0 




0000710 


0062 






00 80 J=1.N 




0000720 


0063 






IF(ZC( J).GE.O.O)GO TO 80 




0000730 


006A 






DO 75 1=1. IM 






0065 






IF( IBASI I ) .EO. J) GO TO 80 






0066 




75 


CONTINUE 






0067 






XE= -WS( J) /ZC( J) 




0000740 


0068 






IF(TEMP.LE.XE)G0 to 80 




0000750 


0069 






TEMP = XE 




0000760 


0070 






IND = J 




0000770 


0071 




80 


CONTINUE 




0000780 


0072 






IF( IND. EO.O) go to 100 




0000790 


0073 






CALL DUAL! IND) 




0000800 


007A 






CALL SIMPLXI IND) 




0000810 


0075 






GO TO 60 




0000820 




L 

C 


CHECK FOR INFEASIABILITY 






0076 


L 


100 


IF(ZCB( 1) )20C,110,120 




0000830 


0077 




110 


IF(ZCB( 2) .LT.0.0) GO TO 200 




0000840 


007R 


c 


120 


IF(IM.GE.M) go TO 190 




0000850 




c 

r 


ADD 


NEXT RESTRAINT 






C079 


L 


130 


IM = IM 4^ 1 




0000860 


0080 






CALL REST 




0000870 


0081 






lUNB = 0 




0000880 


0082 






GO TO 60 




0000890 


0083 




190 


CALL TIMEIT(-l.CL) 




0000900 


0084 






WP I TE( 6.5000) I PROB ,CL .ITER 




0000910 


0085 




SCOOOFHRMAK lOX, 'PROBLEM' , 14, ' TIME IS ' , -(SP F 1 5 . ft , ' 


SECONDS WITH', 


0000920 








116. ' ITERATIONS'.//) 




0000940 


0086 


r 




if(ko.eo.o)GO to 6 




0000990 




c 

r 


K0= 


0 IS FOR USE OF ADDITION OF RESTRAINTS 






0087 






IM=2 




0001000 


0088 






K0 = 0 




0001010 


0089 






GO TO 5 




0001020 


0090 




500 


STOP 




0001030 


0091 




200 


WRITEI6.4010) 




0001040 


0092 






GO TO 190 




0001050 


0093 


4010 


FORMAT! ICX. 'SOLUTION INFEASIBLE'//) 




0001070 


0094 






END 




0001080 



32 



ORTRAN IV G L6VtL 1, MOO 1 



MAIN 



DATE 



6815Q 



05/22/31 



0001 


C 




SUBROUTINE SIMPLX(J) 


0001090 


0002 






00 1 ME NS ION P( 53, 151 ) .B( 52) ,C( 15 It ,BT( 53,2 » , I BAS ( 53 ) .ONE ( 53 ) , TEM( 2 ) , 


0001100 








lPINV(53,53),ZCfl(2),k(15li,WSB(2),WS(151),T(2),X(53l, DBAR (53) 


000111^ 


00C3 






COMMON/ INT/ I UNR, IM,M,N,IBAS,ITFR 


0001120 


C004 






CGMMON/FLOAT/P,B,C,PT,ONE,TEM,PINV,ZCB,ZC,WSP,WS,T,X,OBAR, EP 


0001125 


0005 


r 




lUNB = 0 


0001130 




L 

c 

r 


FOR 


J=1 THE ONLY CHANGE IS THE BASIS AND ZC(J)'S 




00C6 


Lf 

r 


5 


IF(J.EO.l) GO TO 130 


00011 40 




L 

c 


COMPUTE X(JI ANP FIND THETA FOR PRIMAL = M I N ( X ( I , B ) /X < I , J ) ) OVER J 






c 

r 




FOP X(I,J)>0 




0007 


L 




T( 1) = 999999. 


0001 150 


0003 






T( 2) = T( 1 ) 


0001 160 


0009 






DO 20 1=1 , IM 


0001170 


0010 






L = IM-Ul 


0001180 


0011 






X(L) = 0.0 


0001 190 


0012 






00 10 K=1 , IM 


0001200 


0013 






V = PINV(L ,K)^P(K, J) 


0001205 


0014 






X( U = X( L) + Y 


0001210 


0015 




10 


CALL ROUND(X(L) ,Y,FP) 


0001211 


0016 






IF(X(L) .LE.O.)GO TO 20 


000121^^ 


0017 






TEM( 1) = aT(L,n/X(L) 


00012'’0 


0018 






TEM( 2) = BT(L»2> /X(L) 


0001230 


0019 






IF{ T( n-TEM( 1) )2C, 11 , 12 


0001240 


0020 




11 


IF( T( ?) .LE .TfcM( 2) )G0 TO 20 


00012*^0 


0021 




12 


T( n = TEM( 1 ) 


0001260 


002? 






T( 2) = TEM(2) 


0001270 


0023 






ID = L 


0001280 


0024 


r 


20 


CONTINUE 


0001290 




Lr 

C 

r 


COVPUTE NEW TABLEAU 




0025 


L 




IF( T( 1 ) .FQ. 999999. )G0 TO 50 


0001300 


0C26 






IBAS( ID) = J 


0001310 


C027 






ONE( 10) = 0.0 


0001320 


0028 






00 30 1=1 , IM 


0001 330 


0029 






IF( I .EO. in ) GO TO ^0 


0001333 


C03C 






Y = T( 1 H'XI 1 ) 


000133^ 


0031 






BT(I , 1 ) = 8T( I , 1 ) - Y 


0001334 


0032 






CALL ROUNO(RT( 1,1) ,Y,EP) 


0001335 


0033 






Y = T(2)’<‘X( I ) 


0001 336 


0034 






BT( I ,2) = BT( I ,2) - Y 


0001336 


0035 






CALL ROUND! BT( 1,2) ,Y,EP) 


0001337 


0036 






XE = X( I ) /X( ID) 


0001360 


0037 






DO 30 K=1,IM 


0001370 


0038 






Y = PINV( ID,K)’«‘XE 


0001371 


0039 






P INV( I ,K) = PINV( I ,K ) - Y 


0001372 


0040 






CALL R0UN0(PIMV( I ,K) ,Y,EP) 


0001 375 


0041 




30 


CONTINUE 


00013^5 


0042 






BT( 10,1) = PT( in,l)/X(ID) 


0001387 


0043 






BT( ID, 2) = BT( ID,2) /X( ID) 


COOl 398 


0044 






DO 31 1 = 1 , IM 


cool 386 


0045 




31 


PINV(ID,I) = PINV( ID,I )/X( ID) 


0001389 


0046 




70 


ZCB( 1 ) = 0.0 


0001390 


0047 






ZCB( 2) = 0.0 


0001410 


0048 






DO 80 1 = 1 , IM 


0001420 


0049 






ZCB(l) =ZCB(1) ^GNE(I)*RT(I,1) 


0C01430 


005C 




80 


ZCB(2) =ZCB(2) ^-nNF(I)«BT( 1,2) 


0001440 


0051 






ZC( 1) = ONE! 1) 


0001450 


0052 






on 90 L=2,N 


000146C 


0053 






ZC(L) = 0. 


0001470 


G054 






on 90 1 = 1 , IM 


0001490 


0055 






on 90 K=1 , IM 


0001500 


0056 






Y = ONF( I)’«^PINV( I,K)+P(K,L) 


0001505 


0057 






ZC(L) = ZC(L) ^ Y 


0001510 


0058 


c 


90 


CALL ROUNO(ZC<L) ,Y,EP) 


0001515 



33 



FORTRAN IV G LEVEL 1, MOD 1 SIMPLX DATE = 68159 05/22/31 



0059 


e 


CHECK FUR AN ADDITION ITERATION WITH PRIMAL 
TEMP = 0.0 


0001520 


0060 




DO 120 L=1,N 


0001530 


0061 




IFIWS(U.NE.O.O) GO TO 120 


0001540 


0062 




DO no I = 1»IM 


0001550 


0063 




IF( IBAS( I ) .EO.UGO TO 120 


0001560 


006^ 




no CONTINUE 


0001570 


0065 




IF(ZC(L).GE.TEMP) GO TO 120 


0001580 


0066 




J=L 


0001590 


0067 




TEMP = ZC(U 


0001600 


0068 




120 CONTINUE 


0001610 


0069 




ITER = ITER ♦ 1 


0001615 


GOTO 




IF(TEMP.NE.O.O) GO TO 5 


0001620 


0071 




200 RETURN 


0001630 


0072 




130 IBASm = 1 


0001640 


0073 




ONE! 1) = 0.0 


0001650 


0074 




ID = 1 


0001655 


0075 


r 


GO TO 70 


0001660 


0076 


c 

c 


SOLUTION UNBOUNDED 
50 IF( IM.GE.M) GO TO 60 


0001670 


0077 




lUNB = 1 


0001680 


0078 




GO TO 200 


0001690 


0079 




60 WRITE(6,4020) 


0001700 


0080 




lUNB = 2 


0001710 


0081 




4020 FORMAT(10X»*SOLUTION UNBOUNOEDV/) 


0001720 


0082 




GO TO 200 


0001730 


0083 




END 


0001740 
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FORTRAN IV G LFVEL If MOD 1 



MAIN 



DATE = 681S9 



C5/22/31 



0001 

0002 

C003 

OOCA 



000 5 
0006 
0007 

ooca 

OOOQ 

OCIC 

0011 

0012 

0013 

0014 

0015 

0016 
0017 



e 

c 

SUBRDUTINF DUAL(J) 

OD I MENS ION P( 5? , 151 ) , H( 52) , 0115 1) , BT ( 53, 2) , IBAS ( 5^ ) , ONE (51),TFM{2), 
lPINV(53,53)fZCB(2),/C(l5l),WSB(2),WS(151),T(2)fX(53),r)RAP(53) 
COMMJN/INf/I UNB, IM,M,N,IBAS,ITFR 

CCMMON/FLOAT/P, BfC ,BT,0NF,TFM,PINV, ZC B , Z C , WS R , W S , T , X , OB AP , E P 

C 

C CGMPUTF NFW W$(J) WHICH IS EOUIVAl.FNT TH S(J) 

C 

XF = WS(J)/ZC( J) 

DO 10 K=1,N 

Y = ZC( K ) «XF 

WS ( K ) = WS ( K ) - Y 

10 call RnUND(WS(K) ,Y, EP) 

Y = ZCBI 1 ) 

WSB( 1) = WSP( 1 ) - Y 
call ROUND! WSB( 1 ) , Y, EP) 

Y = ZCB(2)*XF 

WSB! 2) = W$B(2) - Y 
CALL PQUNO(WSR(?) ,Y,EP) 

return 

END 



0001750 
0001760 
0001770 
0001 78^^ 
00017B5 



0001790 
0001 aoo 
OOOIROI 

oooiao^ 

0001 BOR 
0001 R04 
0001B05 
C001B06 
0001B07 
0001808 
0001810 
0001840 
0001850 
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FORTR 

0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 



ooo ono ooo 



IV G LEVEL 1, MOO 1 



MAIN 



DATE = 68159 



05/22/31 



e 

c 

SUBROUTINE REST 0001860 

00 1 PENSION P( 53* 151 ) , B( 52) ,C( 151) .BT(5 3,2) , I BAS ( 53 ) * ONE ( 53 ) * TEM( 2) * 000 1870 
1PINV( 5 3,53 ) *ZCB( 2) , ZC( 151 ) *WSR( 2) , WS( 15 1 ) *T ( 2 ) , X ( 53 ) * OB AR ( 53 ) 0001880 

COMMON/INT/tUNB* IM*M,N.IBAS*ITER OOOIRQO 





COMMON/FLOAT/P,B,C,BT,6NE,TEM,PINV,ZC8,2C,WSB,WSf T,X,DBAR,EP 


0001895 


COMPUTE D-BAR AND B VECTOR 






K = IM-1 


0001900 




BT( IM,1» = 0.0 


0001910 




BT(IM,2) = B(KI 


0001920 




DO 10 1=1. K 


0001930 




IF( IBASI I ) .GT.NI GO TO 5 


0001935 




L = IBAS( I ) 


0001940 




DBARI I ) = -P( IM.U 


0001950 




BT(IM.I) = BT(IM.l) DBARI I )*BT( 1 , 1) 


0001960 




BTIIM.2) = BT(IM,2) ♦ DBARI I )*BT(I ,2) 


0001970 




GO TO 10 


0001971 


5 


DBARIII = 0.0 


0001972 


10 


CONTINUE 


0001973 




IFIBTI IM, 1 n 5C.2C.30 


0001980 


20 


IFIBTI IM,2».LT.0.0I GO TO 50 


0001990 


COMPUTE NEW INVERSE AND COST VECTOR 




30 


DO 40 J=1.K 


0002000 




DO 40 1=1, K 


0002010 


40 


PINVIIM.JI = PINVIIM.J) .• DRARII)*PINVII.J) 


0002020 




ZCIII = ONFIl) 


0002021 




DO 45 J=2.N 


000202? 




ZCIJ) =0.0 


0002023 




DO 45 1 = 1, IM 


0002024 




DO 45 IK=I.IM 


0002025 




Y = ONE I I )*PINVI I . IKI*PI IK. J) 


0002026 




ZCIJ) = ZCIJ) ♦ Y 


0002027 


45 


CALL ROUNDIZCI J) .Y.EP) 


0002028 




ZCBI 1) = ZCBIl) - BTIIM.l) 


0002029 




ZCBI2) = ZCBI2) - BTIIM.2) 


0002030 


100 


RETURN 


0002031 


USED TO INSURE FEASIABILITY BY MAKING THE fl COMPONENT NON-NEGATIVE 




50 


BT( IM, 1) = -PT( IM, 1) 


0002040 




BT( IM,2) = -BT( IM,2) 


0002050 




DO 60 1=1, K 


0002060 


60 


DBARI I ) = -DBARI I ) 


0002070 




PINVI IM,IM) = -1.0 


0002130 




GO TO 30 


0002140 




END 


0002150 
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ORTRAN IV G LbVEL 1, MOO 1 



VAIN 



OATF = 681SQ 



05/22/31 



00 Cl 


e 

c 

SUBKOUTINE RCUNOt A,R,EP) 


0002160 




C A=StARTlNG VALUE, R^ADOFf) VALUE, EP= LOWEST ROUNO-nFF OESIPFO 


0002170 


0002 


IF(A.Eg.O.O.rR.R.EO.O.O) GO TO 200 


0002150 


0003 


IF ( ARS( A/B ) .GT .6P I GO TO 200 


0002190 


0004 


A = 0.0 


0002200 


0005 


200 RETURN 


0002210 


0006 


6N0 


00022^0 
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FORTR 

0001 

0002 

0003 

000 ^ 

0005 

0006 

0007 

0008 

0009 

0010 

0011 



ooo 



IV G LEVEL 1, MOO 1 



MAIN 



DATE = 68159 



05/22/31 



E 

C 



SUBROUTINE TIMEIT(N,TIME) 0002230 

N=0 STARTS CLOCK, N=-l STOPS CLOCK 





IT=N+2 


0002270 




GO TO ( 20,10) , IT 


0002280 


10 


CALL TIMONIM) 


0002290 




TIMEM=M 


0002300 




RETURN 


0002310 


20 


CALL TIMOFF(M) 


0002320 




TIME=M 


0002330 




TIME=(TIMEM-TIME)*26.0 


0002340 




RETURN 


0002350 




END 


0002360 
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FORTPAN IV G LEVEL 1, MOD 1 



MA IN 



DATE = 6RI‘59 



OS/22/31 



f 



OOCl 


C 


SUBROUTINE PSEUDO 


C00737C 


0002 




00 I MENS I ON P( S3,151 ) ,B( S2( ,C( 15 1) ,BT( 53, 2) , I R AS ( 5 3 1 , ONE ( 53 1 , TFM( 2) , 


0002380 






lPINV(53,53),7CB(2),/ftl51),WSP(2),WS(l51),T(2),X(53), OBARI 53) 


0002290 


0003 




COMMON /I NT/ I UNR, IM,M,N,IPAS,ITER 


0002400 


COO^ 


r 


COMMON/FLOAT/P.B.C.BT.DNE ,TEM,PINV,?CB, 7C,WSB,WS,T,X,DRAR,rP 


0002410 




L- 

c 


M = NUMBER CF CCNSTRAINTS, N = NUMBER OF VARIABLES INCLUDING M SLACK 






c 


VARIABLES 






c 

r 


maximum M = 50, MAXIMUM N = 150 






c 

r 


CHANGE STATEMENTS 1 AND 2 TO DESIRED PROBLEM SIZF 




00C5 




1 M = 20 




0006 




2 N = 70 




0007 




IM = M <• 1 




0008 




IN = N 1 




OOOo 




K1 = N - M + 1 




OOlC 




K2 = K1 + 1 




0011 




C(J) = -URN(l) 


0002430 


0012 




no 20 J=2,K 




0013 




IFICI Jl .LE.-.OOl ) GO TO 10 




0014 




CIJI = 0.0 


0002450 


0015 




10 DO 20 1=2, IM 




0C16 




P( I , J) = URNI 1 ) 


0002470 


0017 




IFIPI l,J).GT. 0.001) GO TO 20 


0002480 


0018 




P( I , J) = 0.0 


0002490 


0019 




20 CONTINUE 


0002500 


0C2C 




DO 30 1 = 1, M 




00 21 




Bin = URNI 1 )*5.0 


0002520 


0022 




1F(B( II .GT. 0.011 GO TO 30 


0002530 


0023 




Bin = 0.0 


0002540 


0024 




30 CONTINUE 


0002550 


0025 




DO ‘tO J = K2,1N 




0026 




CIJI =0.0 




0027 




DO 40 I = 2,IM 




0028 




40 P( I , J) = O.C 


0002580 


0029 




DO 50 1=2, IM 




0030 




K = N- M ♦ 1 




0031 




50 PI I,KI = -l.C 




0032 




RETURN 


0002640 


0033 




END 


0002650 
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a significant improvement in the solution time by the use 
of adding the constraints one at a time. 
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